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We have studied numerically the evolution and decay of axion strings. These global defects decay 
mainly by axion emission and thus contribute to the cosmological axion energy density. The relative 
importance of this source relative to misalignment production of axions depends on the spectrum. 
Radiation spectra for various string loop configurations are presented. They support the contention 
that the string decay contribution is of the same order of magnitude as the contribution from 
misalignment. 



1. INTRODUCTION 

Axion strings arise in the early universe 
when the global U(1) PQ symmetry breaks 
spontaneously at the scale f a ~ 10 12 GeV. A 
Brownian network of strings is formed initially 
with O(l) string per horizon volume. As the 
horizon expands, wiggles on the strings will 
start oscillating and radiate axions. Strings will 
also intercommute and form closed loops, 
which collapse and convert into axions. 
Simulations of gauge string networks show 
that a scaling solution [1,2] is reached after a 
few Hubble times with about one long string 
per horizon size and a population of decaying 
loops. This process continues until QCD time, 
when the axion acquires a small mass. As a 
consequence, domain walls form between 
strings. The wall-string system is short-lived 
because wall tension pulls the strings together, 
followed by mutual annihilation into free 
axions. If there is an inflationary period after 
PQ symmetry breaking with a reheat 
temperature < f a , strings get diluted and do 
not contribute to the cosmological axion 
density. 

The dynamics of axion strings is governed 
by the classical Lagrangian density 

L=\d^d^-\X{^-flf (1) 

where (j) = |0|exp(ia// a ) is a complex scalar 

field whose phase is the axion field. A simple 
solution of (1) is a static, straight global string 



with core size 8=(y[Xf a ) *and energy per 
unit length 

H « nfZ\n(LlS) (2) 

where L is the long distance cutoff, i.e. the 
interstring distance. The logarithmic factor in 
Equation 2 distinguishes a global from a gauge 
(local) string, whose energy density decays 
exponentially outside of its core. For axion 
strings near QCD time, In (t qcd fa ) ~70 an d 
most of the string energy resides outside of the 
core. 

The number of axions emitted by strings 
during the string epoch (1 GeV < T < 10 12 GeV) 
can in principle be calculated from the 
radiation spectrum of the evolving string 
network. Analytic techniques are not well 
developed and one must rely on numerical 
simulations. A major difficulty is the 
enormous loop to core size ratio of realistic 
axion strings. Currently, the largest computer 
simulations have \n{Ll8)~l. Yet, it is 
valuable to study the spectrum as a function of 
ln(L/<5) in order to deduce a possible trend. 
This is especially important for closed loops as 
they dissipate most of the string network 
energy. 

In the past, there has been considerable 
controversy over the correct radiation 
spectrum. One group [3,4,5,6,7] (scenario A) 
argues that the spectrum is strongly peaked at 
wavelengths of order the loop size. This 
corresponds to an under-damped decay with 



about 5-10 oscillations per loop half-life when 
ln(L/5)=70. A second group [8,9,10] 
(scenario B) argues that the loops decay 
without oscillations, in a time of order the 
initial size divided by the speed of light. The 
radiation spectrum is dEldk~\lk (which is 
also the spectrum of the static string field) with 
cutoffs at the inverse loop size and core 
thickness. Of course, it is possible that neither 
scenario A or B is correct. 

The cosmological axion number density at 
time t is related to the radiation spectrum 
through 



1 f f dt' n{t') 
HaW ~ t 3/2 t ! pe t' 3l2 <o(t') 



where 
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(3) 



(4) 



and dE/dco is the radiated axion energy 
spectrum. A useful quantity is 



configurations are set up on large (~10 7 points) 
Cartesian grids, and then time-evolved using 
the finite-difference equations derived from 
Eq. (1). A FFT spectrum analysis of the kinetic 
and gradient energies during the collapse 
yields N ax (f) . 

2.1. Circular loops 

Because of azimuthal symmetry, circular loops 
can be studied in r-z space. By mirror 
symmetry, the problem can be further reduced 
to one quarter-plane. 

The static axion field far from the string core 
is (with f a - 1) 



a (r,z) ~ —Q.(r, z) 



(6) 



in the infinite volume limit, where £1 is the 
solid angle subtended by the loop. We use as 
initial configuration the outcome of a 
relaxation routine starting with Equation 6 
outside the core and 



„ . . ,dE . .dk 

Nax(t) = \ — (t) — 

dk k 



(5) 



0(p) = O.58(p/S)exp(/0) 



(7) 



In scenario B, N ax is predicted to stay constant 
since both axion radiation and static string 
field have the same spectrum. Conversely, one 
expects N ax to increase by ln(LA5) in scenario 

A. The factor r = A^ nal / A?™ tial , by which N ax 
increases during the decay of a string loop 
determines the string decay contribution to the 
axion cosmological energy density. The string 
decay contribution is r times the contribution 
from misalignment. In scenario A, r is of the 
order ln(L/c5) , whereas in scenario B, r is of 
order unity. 



2. LOOP SIMULATIONS 

We have performed simulations of various 
loop geometries: (1) circular loops initially at 
rest, (2) noncircular loops with angular 
momentum, and (3) string-antistring pairs 
with angular momentum. The initial 
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Figure 1: Core position of collapsing loop versus time for 
X = 0.001 (dotted line) and X = 0.004 (solid line). The 
lattice size is 4096x4096. 

within the core. Here, p is the distance to the 
string, and 6 is the winding angle. The 
relaxation and the subsequent dynamical 
evolution are done with reflective boundary 



conditions. A step size dt = 0.2 was used for 
the time evolution and the total energy was 
conserved to better than 1 %. In general the 
loops collapsed at nearly the speed of light 
without a rebound. For a small range of 
parameters, 80 < R /8 < 190, where R is the 
initial loop radius, we noticed a small bounce 
as shown in Figure 1. 

There is a substantial Lorentz contraction of 
the string core as it collapses (see Figure 2). A 
lattice effect became evident when the reduced 
core size becomes comparable to the lattice 
spacing. This lattice effect consists of a 
"scraping" of the string core on the underlying 
grid, during which the kinetic energy of the 
string gets dissipated into high frequency 
axion radiation. We always choose X small 
enough to avoid this phenomenon. 
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Figure 2: Intensity plot of potential energy (contours 
represent constant potential energy) in vicinity of string 
core for R,= 2400, X = 0.001 at f = 2600. The Lorentz factor is 
about 4. The arrows represent the axion field. 

A spectrum analysis of the fields was 
performed by expanding the gradient and 
kinetic energies as 



v z = Lcmn Jo(k m r) sin(k n z) 

mn 

V r <t>= LCmn J i (k m r) cos(k n z) (8) 

mn 

<P= Y.c mn J Q (k m r)cos(k n z) 

mn 

with the boundary conditions /j (k m R max ) = 
and sin (k n Z max )= . The dispersion 

/ 9 9 

relationship is given by (0 mn = -^k m +k n . 

Figure 3 shows the power spectrum 
dE/d(lnk) displayed in Ink bins of 
width A In k = 0.5 at t = and after the collapse 
at t = 3000. At both times, the spectrum 
exhibits an almost flat plateau, consistent with 
a dEI dk^Xlk spectrum. The high frequency 
cutoff of the spectrum is increased however 
after the collapse and is associated with the 
Lorentz contraction of the core. 
The evolution of N ax = T.(E mn lk mn ) during 

mn 

the loop collapse was studied for various 
values of Rq/8 . We observe a marked 
decrease of by ~ 20 % during the collapse 
roughly independent of R Q 15 . 
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Figure 3: Energy spectrum of collapsing loop for R„ = 
2400 and X = 0.001. The white (black) histogram represents 
the spectrum at t = (f = 3000). The increased high 
frequency cutoff of the final spectrum is due to the Lorentz 
contraction of the core. 
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Figure 4: N ax of a circular loop as a function of time for 

R,= 2400, and X = 0.004 (solid line), X = 0.001 (dashed line), 
and X = 0.00025 (dotted line). 

2.2. Rotating loops 

There exists a family of nonintersecting 
("Kibble-Turok") [11,12,13,14] loops, which 
have been studied in the context of gauge 
strings. They are solutions of the Nambu-Goto 
equations of motion and all have non-zero 
angular momentum. 

Intercommuting (self-intersection with 
reconnection) causes the loop sizes to shrink, 
and hence the average energy of radiated 
axions to increase and hence to decrease. 
Intercommuting favors scenario B for these 
reasons. We picked the Kibble-Turok 
configuration as an initial condition to avoid 
intercommuting as much as possible, thus 
giving scenario A the best possible chance to 
get realized. 

A common loop parameterization is given by 

x(a,t) = -|{(l-a)sin<7_ + -^asin3c7_ +sinc7 + } 

y(a,t) = y{-(l-ai)cosc7_ --i-acos3a_ -cosi/f cos<7 + } 

z(o,t) = ^ {- 2 ^/a(l - a) cosct. - sini//xosc7 + } 

(9) 

where a+ = {a ± t)IR , and ae (Q,2nR) is the 
length along the loop. For a significant subset 
[12,13,14] of the free parameters cue (0,1), 
\jfB{—K,n) the loop never self-intersects. A 
noteworthy feature is the periodic appearance 
of cusps, where the string velocity 
momentarily reaches the speed of light. The 
motion of a gauge string is damped by 



emission of gravitational radiation and the 
loop diameter shrinks with time. The power is 

P = TGji , where /J. is the energy per unit 
length, G is the gravitational constant, and r 
is a constant ~ 50 - 100 which depends on 
a, y/ . Evidently, the power is independent of 
loop size, and the loop undergoes 
l/(TG/J.) =10 4 oscillations in its lifetime. The 
radiation power spectrum was numerically 
determined [15,16] to be P k ^ l/k 4n . 

Axion strings are much more short-lived 
than gauge strings and radiate axions 
efficiently due to the strong topological 
coupling between string and field. No closed 
loop solutions are known however. 

According to scenario A [7,15], the radiated 

power is P = 50 f a , independent of loop size, 
and the loop shrinks linearly with time. The 
expected number of oscillations per loop half- 
life is 7Tln(L/<5)/50 = 4 for In (L/S) =* 70 . In 
addition, N ax should increase by a factor 
In {LI 8) . In scenario B on the other hand, N ax 
should remain approximately constant. 

We performed numerous simulations of 
rotating loops on a 3D (256 3 ) lattice with 
periodic boundary conditions. Standard 
Fourier techniques were used for the spectrum 
analysis, and N ax was computed as a function 
of time using the dispersion relationship 

0) mnp = ^2(3 - cos k m - cos k n — cos k p ) . (10) 
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Figure 5: N at of non-intersecting ("Kibble-Turok") loops 

as a function of time for a = 0.01, <|) = 0, and X = 0.2 (solid 
line), X = 0.1 (dashed line), and X = 0.0625 (dotted line). The 
lattice size is 256 3 and R = 72. 
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Figure 6: Energy of non-intersecting loop. Shown are 
total, gradient, kinetic, and potential energy as a function of 
time for a = 0.01, \|/ = 0, and X = 0.0625. 

Figure 5 shows N ax for various core sizes 
and constant a, y/ . The behavior is very 
similar to that of a non-rotating circular loop 
with a reduction of ~ 25 %. Figure 6 depicts the 
energy of the collapsing loop. Clearly, the total 
energy is well conserved. A few percent of the 
loop energy is dissipated as massive radiation, 
shown here as E pot . 

2.3. String-antistring pairs 

Lastly, we studied rotating, parallel string- 
antistring pairs, which can be thought of as 
cross-sections of loops with large eccentricity. 
The 2048 2 lattice was initialized with the 
Abrikosov ansatz [17] 

4>{x,y) = 4> l {x-x h y-y l )4>l(x-X2,y-y 2 ) (11) 

where (x i ,y i ),(x 2 ,y2) are the locations of 
string and antistring respectively. The fields 
were relaxed with periodic boundary 
conditions and the cores held fixed. The time 
derivative <j)(x,y) was obtained by a finite 
difference over a small time step. The 
attractive force per unit length between strings 



system to form, which decays by emission of 
axions and eventually annihilates. 
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(12) 



is competing against the centrifugal force. In 
general, one expects a bound relativistic 



Figure 7: Collapsing string-antistring pair. The initial 
rotational velocity of the cores is 0.4 and X = 0.005. The 
snapshots are for t = 40, 800, 1440. 



Figure 7 shows the decay of a string-antistring 
system for an initial rotational velocity of 0.4. 
For the parameter range (3 < ln(p/(5), 
v rot < 0.6, p ~ 500), the system collapsed 
immediately in a time of order p . 

The spectrum and N ax were computed as 
well and Figure 8 shows the evolution for 
several cases. 
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Figure 8: N ax of collapsing string-antistring pair. The 

initial rotational velocity is 0.6 (dashed), 0.4 (solid), 0.2 
(dotted) and X = 0.005. 

3. CONCLUSION 

Two different opinions exist about how 
axion strings decay. The first view (scenario 
A) states that the loops are moderately 
damped and oscillate several times during 
their lifetime. The instantaneous radiation 
spectrum is dominated by wavelengths of 
order the loop size, and N ax is predicted to 
increase during the decay by a factor In {LIS ) . 

Scenario B on the other hand proposes that 
the loops are critically damped, the radiation 
spectrum is 1 / k , and N ax is approximately 
constant. Our simulations clearly support 
scenario B since N ax decreases by about 20 % 
during the decay of axion loops for In {LIS) ~ 
6. Moreover, the decrease of N ax which we 
observe during the decay of string loops, does 
not change with \n(Ll8) over the range of 

In {LIS) which we investigated. 
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